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Abstract 

This paper settles the exister\ce question for a rather general class of convex optimal design 
problems with a volume constraint. In low dimensions, we prove the existence of an optimal 
configuration for general convex minimization problems ruled by bounded measurable degen- 
erate elliptic operators. Under a mild continuity assumption on the medium, the free boundary 
is proven to enjoy the appropriate weak geometry and we establish the existence of an optimal 
design for general convex optimal design problems with volvime constraints for all dimensions. 
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1 Introduction 

Well known for modeling important problems in applied mathematics, respected for the chal- 
lenging mathematical questions they give rise to and admired for their intrinsic beauty, opti- 
mization problems with volume constraints have received an overwhelming attention in the past 
few decades. In general, the usual techniques of the Calculus of Variations are not sufficiently 
powerful, or even appropriate, to establish existence of optimal configurations for those classes 
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of problems. This fact has inspired remarkable recent advances in a number of branches of ap- 
plied analysis in an attempt to develop the right set of analytical and geometrical tools to study 
optimal design problems with volume constraints. 

One of the fundamental motivations of this present work can be, in its most basic form, stated 
as follows: given an 71-dimensional body and a fixed amount of insulating material, what is the 
best way of insulating it? Depending on the flexibility allowed, the mathematical set-up used to 
model this classic question can also be employed in the analysis of a variety of other problems 
in applied mathematics. In more precise mathematical terms, but still using the language of heat 
conduction, the above question takes the following form: let D he a fixed Lipschitz boimded 
domain in R" (the body to be insulated), (p: dD — > M be a prescribed positive function (the 
temperature distribution on D), and t > be a given positive number (the amount of insulating 
material available). For each configuration Vl that surrounds D and obeys D) < l, we 

compute the flux associated to it: 

In general, 3 is related to a boundary integral involving a potential uq, linked to by a prescribed 
PDE. The optimal design problem is then 

Min {Z{n) \ ^^D and £"(fi \D)<i] . (1.1) 

Probably the first and still one of the most influential works in this line of research is the pio- 
neering article of Aguilera, Alt and Caffarelli, IAAC86I . In this paper, the authors address the 
question of minimizing the Dirichlet integral when prescribed the volume of the zero set. More 
precisely, they study the optimization problem 

Min iVupdX I u e ffi(f7), u = Lp>Qondn and ({u = 0}) = a| , (1.2) 

for a fixed a < iL"(r2). In the case of an exterior domain, O = M" \ D, problem (|1.2^ can be 
used to model a very simple, yet interesting optimal design problem with volume constraint as 
stated above. Namely, suppose D is evenly heated. If one tries to minimize the heat flux given by 
Idn ""/j*^^""^ (^)' where u is the capacity potential associated to 51, with iL" {il \ D) prescribed, a 
simple application of Green's identity reveals that the heat flux equals the Dirichlet integral, and 
therefore the problem becomes identical to lll.2b . Fine regularity properties of the free boundary, 
d{u* > 0}nil, where u* is a minimizer of l|1.2l l rely on the powerful geometric-measure machinery 
developed by Alt and Caffarelli in IAC81I : the magnum opus of free boundary regularity theory 
for variational problems. 

A significant generalization of problem (|1.2| | was carried out by Lederman in lLed96l . In 
this paper, the author studies the non-homogeneous minimization problem, that is, the Dirichlet 
integral is replaced by / \S/u\'^dX — / gu, for a given g bounded away from zero. 

In an important paper, Ambrosio, Fonseca, Marcellini and Tartar, [AFMT99], address another 
major generalization of problem \1.2\ . Namely they establish the existence of a minimizer to the 
functional J :— W{Du)dx, for W : M''^" (0, 00) and quasi-convex, with the multiple 
volume constraint = Zi}) ~ ai,0 < i < k. In a subsequence article, Tilli, in |Tilli99J, 

showed, for W{^) := |^|^, that in the case of just two level constraints, the minimizers are locally 
Lipschitz continuous. 

Still assuming a constant temperature distribution, Oliveira and the author in IOT06I studied 
the optimization problem jl.lh governed by the p-Laplacian operator when the flux is given by 
J(u) := Jgj2 {u^Y~^ dTC"~^{X). This translates into the analysis of the minimization problem 
l|1.2b , for the p-Dirichlet integral, that is, W{£_) ^ for p > I. 

The first work to deal with optimal design problems with non-constant temperature distri- 
bution (p: dD (0,00) is IIACS87I . In this paper, the authors consider the linear functional: 
3(f2) = J AudX, where u is the harmonic function in fl \ D, taking boundary data ip on dD 
and zero on 90. Even for this simple fimctional, major difficulties arise. For instance, the free 
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boundary condition, that is, the behavior of Vu* along the free boundary, dfl*, is non-local and it 
required a new machinery to establish the appropriate geometric-measure properties of the free 
boundary necessary to perform suitable smooth perturbations. The latter is used in its entirely to 
finally conclude the existence of an optimal design. 

At least for smooth competing configurations, 17, for the linear functional studied in IIACS87I 
we have 

a(f7) := / AudX = / Uyd7e-^{X)= j u^,d7e-^{X). 
J Jdn JdD 

This is a naive, yet important observation, as the latter integral is taken over the fixed boundary. 
Therefore, at least in an intuitive perspective, a non-linear theory for this class of minimization 
problems should use /^^ u^dH^^^^{X) as its linear pattern. From the applied viewpoint, if one 
allows a nonlinear flux, Z that might also depend upon the local structure of the boimdary of the 
body D, i.e., 

m)-^ [ r{X,u^{X))dK'-^{X) (1.3) 

JdD 

the mathematical model l|l.lb would address several other physical situations, such as: optimal 
configurations in electrostatics, problems in material science, flux dynamics, among many others. 
This nonlinear setting, however still only for problems governed by the Laplacian operator, has 
been studied by the author in BTeixOSI and IITeix07l . 

In this present paper, we settle the existence theory for optimal design problem l|l.Hl with 
nonlinear functionals as in (|1.3l l, when un is linked with 51 by a rather general class of degenerate 
elliptic PDEs. In terms of applications, it greatly extends the range of physical systems that can 
be modeled by this set-up. From the mathematical viewpoint, this project brings a number of 
new rather challenging difficulties in its analysis and modern solutions to various issues com- 
monly found in free boundary problems are developed throughout the paper. Free boundary 
regularity theory for uniform elliptic operators in divergence form with merely Holder continu- 
ous coefficients is currently being developed in order to establish C'^''* smoothness of an optimal 
configuration, up to a possible negligible singular set, |Teix-Prep) . 

The article is organized as follows: in section |2l we describe all the mathematical elements 
involved in the model and the optimization problem is accurately stated in that section. Still in 
section|2l we introduce weak formulations of the optimal design problem that are somewhat 
simpler to be tackled from the mathematical perspective. Basic properties of the functional to be 
minimized are established in section |3l The first existence theorem for a weak formulation of the 
original optimization problem is delivered in section ID In section |5j by letting the penalty term 
blow-up, we establish the existence of an optimal configuration to the optimal design problem 
with volume constraint (|l.lb ruled by totally discontinuous degenerate elliptic operators. For that 
though, a technical restriction on the dimension is necessary. In section|6j under C"^ regularity on 
the medium, a series of results concerning the weak geometric properties of the boundary of an 
optimal configuration to the weak formulation of the original problem i ll. Il l are achieved. These 
are used in section[7]to ultimately derive existence of an optimal configuration in all dimensions. 

2 Mathematical set-up 

Throughout the paper, D denotes a fixed Lipschitz boimded domain in R", ip: dD ^ R is a 
prescribed positive function and i > is a given positive number. Our medium deformation will 
be expressed by A : D'^ x R" R", a measurable p-degenerate elliptic map, that is, 

(a) for each ^ G R", the mapping X ^ A{X, ^) is measurable. 

(b) For a.e. X e D'^', the mapping ^ A{X, £_) is continuous. 

(c) There exists constants < A < A < oo and ap > 1, such that, for a.e. X e D*-^ and all ^ G R", 
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(i) ^(X,e).e> A|er 

(ii) < A|er^ 

(iii) {A{X, ei) - A{X, 6), - 6) > 0, whenever a 7^ 6 and 

(iv) A{X,aO=a\a\P~^A{X,£). 

A typical example to keep in mind is 

A{X,0 = A{xmP-^t 

with ^ bounded measurable, which gives rise to the theory of optimal shape problems governed 
by the p-Laplacian in a totally discontinuous medium. 

Our optimization problem is then formulated as follows: for each domain Vl <Z D satisfying 

{n\D)< L, (2.1) 

we consider the ^-potential, u = u{VL), with the prescribed boundary value (p on the fixed bound- 
ary dD, associated to il, i.e. the unique solution to 

l.u:=AiY{A{X,Du)) = QwiVL\D 

u = ipondD (2.2) 
u = on d^l 

and compute 

Z{n) := / r {X, dAu{X)) drC-^{X) (the flux: quantity to be minimized). 

JdD 

Here r : dD x M — > M is a given function, whose properties will be described soon, and 

dAu{X) := {A{X, Vu{X)),^i{X)) (2.3) 

where ji denotes the inward normal vector defined a.e. on dD. The optimal design problem 
we are interested in is the following: 

Minimize | 17 D Z) and {n\D) < t|. (2.4) 

The analytical (and naturally mild) properties assumed on the nonlinearity V are: 

1. For each X e dD fixed, T{X, ■) is convex and increasing. 

2. For each f e M fixed, dtT{-, t) is continuous. 

T(Y t] 

3. If r(Xo, to) then V{Y, to) VF G dD; otherwise ^ ' \ < L, for a universal constant 

T(X, t) 

L > 0. 

Notice that from 1 the following coercivity condition holds: 

lim / T{X,t)drr-^{X)^+oo. (2.5) 

If 1/) is a positive continuous function defined on dD and 7 is a increasing convex function, 
then 

T{X,t) = i^[X)^{t) 

gives a typical nonlinearity that fulfils the above properties. As in the Calculus of Variations, T is 
chosen based upon the particular problem we are trying to model and no relation whatsoever is 
imposed upon the nonlinearity V and A. 
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Sometimes it is convenient to use the language of heat conduction theory to describe the 
elements involved in our analysis. Thus, D is the body to be insulated, ip represents the tempera- 
ture distribution on dD, l corresponds to the maximum amount of insulating material available, 
5 plays the role of the (generalized) heat flux, which is the quantity to be minimized, and A de- 
termines the inhomogeneous and complexity features of the medium. However it is important 
to highlight that this model is widely applicable to several other situations beyond the bounds of 
the classical heat conduction theory and other interpretation of the model might provide different 
insights on what is reasonable to expect to hold. 

It is noteworthy to point out that, since we are not forcing any regularity assumption on the 
medium A, in principle just Holder continuity estimates are available for an ^-potential u = 
u{Q,). Thus, the ^-normal derivative of u, d^u, as entitled in l|2.3|l is not properly defined. Some 
of our primary results concerning geometric properties of the free boundary will not depend 
upon any smoothness condition on the medium. However, just to grapple with this technical 
inconsistence, we will assume throughout the paper that there exists a small 1 >^ (5o > 0, such 
that 

(i) A is Holder (or even only Dini) continuous in Ds^ ■= {X e R" | dist(X, dD) < Sq^ and 

(ii) ip : dD ^ M is accordiantly smooth. 

(2.6) 

Once more we emphasize that for the first part of this project, condition (|2.6b plays merely a 
technical role and, for sake of applications, it should not be seen as a constraint. 

2.1 Penalty Method and weak formulation 

From the mathematical point of view, the minimization problem l l2.4t carries too many difficulties 
to be approached directly. Instead, we will employ a fruitful penalty method in order to formulate 
weak versions of problem l|2.41 . Such a technique has been successfully employed to study a 
variety problems in applied mathematics. 

The intuitive idea behind a penalization strategy is the following: suppose our problem has an 
"undesired" (from the mathematical perspective) constraint on the competing configurations (in 
our case a volume constraint). We then allow any configuration to compete; however we "charge 
a fee" for those configuration that do not obey the previously set constraint. We expect that, if 
the fee is too high, optimal configurations will indeed prefer to satisfy the original constraint. 

Still in a philosophical perspective, one should expect that an optimal configuration, fl*, of 
problem I I2.4I I satisfies 

{n* \D) = L. 

For that, think of l as the budget available and il* as the ultimate object to be built up. Mathe- 
matically, this fact is indeed easily justified. For instance suppose, for an optimal configuration 
n*, we had 

{n*\D) <L- e, 

for some e > 0. Let Xq e dil* be a free boundary point and p > so that ujnp" < £■ Consider 

n := n*UBp{Xo). 

Thus, competes with il* in the minimization problem l|2.4b and, because of maximum principle, 
u(r2) > u{n*). Taking into accoimt that F is increasing and applying Hopf maximum principle 
on dD, we would conclude 

which contradicts the minimality property of O*. Our conclusion is that in problem ||2.4|| we can 
regard the condition £" (f7 \ D) < l as L" {il \ D) = l. For future reference, let us state this as a 
Lemma. 
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Lemma 2.1. Let Q* be a minimizer of problem (12. 4P . Then £"(il* \D) = l. 

Another general comment: we will always extend the ^-potential u{Vi) by zero outside fi. 
Thus, in the distributional sense, 

L [u{n)] = 0, in f7 = {u{n) > 0} and L [u{n)] > 0, in K" \ D. (2.7) 

Returning to the penalty technique issue: we shall borrow the simple, yet quite clever penalty 
term suggested in ||Tilli99l , that is, for each A > 0, we will consider the penalization term g\ : R+ — s 
IR+, defined by 

gx{t):=X{t-L)+. (2.8) 
We then define the A-perturbed functional, Ja, to be 

f r{x,dAu{x))dn''-\x) + gx{ti''{n\D)). (2.9) 

JdD 

Once more, the idea is the following: we allow Z\ to act on any configuration il D D and, when A 
is big enough, we hope that an optimal design il\ for Z\ will satisfy (il^ \D) ~ l, thus it will 
also be a minimizer for our original optimization problem with volume constraint. Our initial 
goal is then study existence and geometric properties of the penalized problem: 

(*Pa) Minimize among all sets f7 D Zjj. (2.10) 

However, even the penalty problem | |2.10l l is, in principle, too hard to be directly approached. 
Thus, for the time being, it will be more appropriate to initially deal with a weak formulation of 
problem | |2.10|| , which we start describing now. Let (5o be the technical number in l|2.6|l . For each 
6 <K So, we us define the functional set 

V{6) := {/ e W^'P{D^) I / = on / > 0, £/ > 0, £/ = in D5} . (2.11) 

Then we define the sample functional set: 

V U V{6) (2.12) 

<5\0 

and the weak formulation of problem | |2.10l l can then be stated as 

CPr^k) ^^-^i f r{X,dA.fiX))dH''-\X) + gx{L-{{f>0}\D))]. (2.13) 

3 Basic functional and analytic properties 

In this section we establish all the basic and necessary properties on the mathematical elements 
of the problems we are interested in, namely, problems l|2.4b , | |2.10| | and | |2.13l l; however we will 
mostly be concerned with the latter, as it the the weakest formulation among them all. 

We start by stating, as a lemma, a simple yet crucial observation regarding the measure theory 
involved on our optimization problems. The proof is somewhat long, but rather standard and 
we omit it here. 

Lemma 3.1. Let f £ V, as in il.lZi . Then £/ defines a nonnegative Radon measure, in . In 
particidar, for any tp e C{D'^) n W^'P{D^), 

[ ij{X)dfif{X)+ f {A{X,Df),DiiX))dX= I tlj{S)-dAf{S)dH"-\S), 

JB,.(Y) JBr(Y) JdBr{Y) 
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for almost allO<r < dist(y, dD). Also, iftjj G Wq''{D^), there holds 



Furthermore, 



{A{X,Df),D^:{X))dX= / i,{X)d^if{X). 

DC J DC 



lif{W^\D)^ / dAf{X)dH^-\X). 

JdD 



(3.1) 



(3.2) 



Another useful results that we make use throughout the paper is (for a proof in the case of the 
p-Laplacian we refer, for instance, to IIDP05I , page 100): 

Lemma 3.2. Let O be a domain in R" and f G W^'P{0). There exists a constant c = c{n, A) > 0, such 
that 



I \\/{,f -h)\PdX 

{{A{X, Df),Df) - {A{X, Dh),Dh)) dX > c \ 

«(/) 



Vif-h)\PdX 



2/p 



if p>2 

if l<p<2. 



where 



«(/) := 



1 1- 



\Vf\PdX 



and h is the A-harmonic function in O that agrees with / on dO. 

Our first Proposition provides an energy estimate for a minimizing sequences to our opti- 
mization problems. More precisely, we have: 

Proposition 3.3. Let Uj be a minimizing sequence for the functional ^ \. Then, 

I|Vuj||lp(£,c) < C, 
where C depends only on dimension. A, D, ip and V. 

Proof. Let [) = ()p be the p-harmonic function in D'~^ that agrees with on dD^ , that is the solution 
to 

Api) = OinD'=' 

I) = ifondD^ (3.3) 
[, e iyi.p(£)C). 

From the maximum principle, there holds 

< t) < sup Lp. 

dD 

For sake of notation convenience, let us denote LujdX := = Mj/ in Lemma [3. II We clearly 
have 



(3.4) 



(f) - = J {A{X, Duj), ([) - Uj))dX 

= j {A{X, Du), Dh{X))dX - j {A{X, Du), Du)dX. 
From the degenerate ellipticity of A, we can deduce from | |3.4| | that 

\ j \Duj{X)\PdX < j {\)-u)d^ij - j {A{X,Du),D\))dX 

< sup(^-/ij(R" \D) + A / \Duj\P''^\D[)\dX 

dD J 

< sup(p-/ij(R" + ^ / \DujY'dX + Ci I \Di]Y'dX. 

dD 2 J J 



(3.5) 



In the last step we have used Young's inequahty and Ci ~ e^^/p where e satisfies eP/(P~i) = 
pX/2{p — 1). In view of i3.2} and the estimate in l|3.5b , we reach the conclusion that there exists a 
constant Ci, depending only on A, D and (p, such that 

\Vu,\\l^^^a) < Ci (^1 + ^ j^Jj,u,{X)<m''-\X)^ , (3.6) 



where a ^{dD). From the monotonicity and convexity properties of the non-linearity F, 

we derive, for each Y G dD fixed, that 

2F (y, ||v^i,||^,(^e)) <c, + r(Y,^ dAU,{x)dn--\x)^ , 

where C2 is a constant depending only on A, D, tp and F. Once more using the convexity of 
F(y, •), it follows from Jensen's inequality that 

2t(y,\\^u,\\1^ c.)<C^ + - I T{Y,dAU,{X))dH''-\X). (3.7) 

Integrate inequality l|3.7|l with respect to Y over dD and taking into account property (iii) of the 
non-linearity F, we derive 

^^F(r,||Vw,||^,(^c))rfH"-i(y) < C3 (^i + ^^F(x,a^u,(x))dH"-i(x)^ , (3.8) 

where again C3 depends only upon A, D, ip and F. Finally, (13.81 1 and the coercivity of the function 

[ F(X,t)dH""^(X), 

JdD 

see (12.5b , together complete the proof of the Proposition. □ 

In view of the energy estimate provided in Proposition |3]3l it becomes natural to investigate 
the behavior of Z\ over weakly convergent sequences in VK^-^. In this direction we have 

Lemma 3.4. Let fj e H^{D^) be a sequence of functions satisfying, Lfj > in the distributional sense 
and, for some S > 0, Lfj = in Dg :— {X e D'~^ | dist(X, dD) < S] . Assume fj converges weakly to 
f in W^'^{D'^'). Then, Lf >0 in the distributional sense, Lf = in Ds and furthermore 

m + gx (|{/ > 0}|) < liminf Uf,) + (|{/, > 0}|) }. 

Proof. The fact that Lf > in the distributional sense follows easily. In fact, for any nonnegative 

^ e Cg (Z?'^), we have 

(L/, V) - / {AiX,Df),Dij)dX = - lim / {A{X, Df,), D^)dX > 0, 

since, for any j, lim J {ADfj, Dip)dX < 0. A similar computation shows that L/ = in Ds- 

Let us turn our attention to the T4^^'p'-weak lower semicontinuity of the functional Firstly, 
the volume penalty term of the functional ^\ is indeed weak lower semicontinuous, since, up to 
a subsequence, fj{X) f{X) for a.e. X G D"-^ . Thus, by Fatou's Lemma 

|{/ >0}|< liminf |{/,>0}|. 

Since, the penalty factor g\ is non-decreasing and continuous, there holds 

PA (|{/>0}|)< liminf (|{/, >0}|), 
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as desired. We now focus our attention on the functional 3 (w) = /^^ r(X, 9^t;)c?H"~^(X). As 
in the Calculus of Variations, in order to establish the W^'^-weak lower semicontinuity of 3, we 
shall explore the convexity assumption on r{X, ■). Indeed, we start by analyzing functional with 
piecewise linear potential, i.e., functionals with this particular profile: 

dm{v):= [ F^{X,dAv)d7e~\X), (3.9) 

JdD 

where F,„ is of the form 

i^„(X,t)= max {Bfc(X)t + Cfc(X)}, Bk,CkeC{dD). (3.10) 

l<A;<m 

We then label, for each A: = 1, 2, • • • , m, the sets 

Vk{f) := {X^dD\ F„,{X,dAf{X)) = Bk{X)dAf{X) + Ck{X)}. 

rn 

Thus dD = [j T>k{f), and we may assume that Vkif) n = 0, whenever k ^ i. Also, recall 

k=l 

that L/j and L/ define Radon measures in , and since fj f in W^'^, by standard elliptic 
estimates, we have 

t./, ^ L/, 

in the sense of Radon measures. Therefore, using a representation as in i3.2\ , we obtain that, for 
any continuous function ( G C{dD), 



aX)dAl{X)dK'-'{X) < liminf / CiX)dAfj{X)dn''-\X). 
dD JdD 

With the above at hands, we estimate 

m „ 

Srnif) = V/ {Bk{X)dAf + CkiX)}dn''-\X) 



< liminfV / {Bk{X)dAfj+Ck{X)}dn"'\X) 



< liminf 

In other words, we have proven functionals as in | |3.9l l are VK^ ^'-weak lower semicontinuous. 
Finally, under the assumption that T{X, •) is convex we know that for each X G dD there exits a 
sequence of functions Fm{X, t) as in ll3.10l l such that, for any t, 

T{X,t)^ lim F^{X,t). (3.11) 

rn — *oo 

As a combination of (|3.111 and the VF^-^'-weak lower semicontiniiity of each the Lemma 
follows. □ 

The results proven in Proposition 13.31 and in Lemma 13.41 are important piece of information 
towards establishing the existence of an optimal shape for problem ||2.10| |; however, at this precise 
stage, those are not enough. We would like to invite the readers to make a small pause in order 
to appreciate the intrinsic difficulty involved in proving the existence of a minimal configuration 
to the penalized problem l|2.10|l . 

Following the natural scheme, one considers a minimizing sequence, ftj, to the fimctional ^x, 

i.e., 

■dx(ilj) — > mm.jA- 

■' flDD 
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If Uj denotes the ^-potential associated to the configuration flj, it follows from Proposition 
13.31 that, up to a subsequence, Uj converges weakly and almost everjrwhere to a function u G 
W^'P{D'~') which is non-negative. As a consequence of Lemma |334l we have that hu > 0. Further- 
more, 

/ T{X,dAu)dn"-\X) + gx{\{u > 0}\) < minZx. 

JdD ^^^^ 

Therefore, a natural candidate for an optimal shape to problem | |2.10| | is 

17 := {X e ]R\D I u{X) > 0} . 

However, with the information we have so far, it is not possible to guarantee that {^,u) is an 
admissible pair, i.e., that u is the ^-potential associated to fl, or equivalently that 

Lw = in Q. 

In fact, it is not true, in general, that if an ordinary sequence of functions Uj, satisfying Auj ~ 
in {uj > 0}, converges weakly in to u, then Au = in {m > 0}. As a general comment, the 
above described difficulty is one of the features that makes problems with varying domains (free 
boundary problems) notably more deUcate. 

4 Existence of minimizer to problem (^a) 

Well, the scheme presented at the end of the previous section is not pointless: since our sequence 
is converging to a special configuration, namely a minimizer for the functional Ja, we should 
keep the hope that this strong additional ingredient will assure that in fact Lm = in {u > 0}. In 
this section, we will carry this delicate analysis out, which will ultimately allow us to conclude 
problem (|2.13l l has always a minimizer. As we will see, even the weak formulation of the penalty 
version of our primary goal presents rather delicate mathematical issues. This is due in part to 
the adverse environment generated by the non-linear and degeneracy features of A, and in part 
to the non-local structure of the problem. The latter makes local perturbations inefficient, and 
thus more creativity is needed to furnish appropriate competing configurations. 

As for our first result towards the existence of a minimizer for problem (|2.13t , we will pro- 
vide an a-priori estimate on the distance from the free boundary to the fixed boundary. This 
is an important supporting result as it allows to seek for minimizers in a more suitable class of 
configurations. 

However, in order to accomplish such a result, we initially need to study an auxiliary free 
boundary problem in the spirt of IAC81I , which we present now. 

Theorem 4.1. Let Obea domain in R" and O ^ M a nonnegative function. Let Abe a p-degenerate 
elliptic map and assume A{-,0) £ C"^ for all O. Then, for any constant r > 0, there exists a minimizer 
V = Vt to the problem 



Minimize 




Furthermore, v is nonnegative, Lipschitz continuous and nondegenerate away from the free boundary 
d{vr > 0}. 

With the free boundary technology available nowadays, it is not hard to establish the existence 
as well as optimal regularity and nondegeneracy of a minimizer to the above problem. Basically 
there are two procedures that lead to these results: one can directly approach the minimization 
problem, by mixing the strategy as in IIAC81I and IDP05I . Another charming and fruitful strategy 
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is to employ a regularizing technique method, basically by mixing the estimates in IIMT07I or 
ITeixll and | DPS0 3j, see also IIK06I . Mathematically the latter is described as follows: choose your 
favorite nonnegative boimded real function (3, such that supp/3 = [0, 1] and, say, l3{Q)dC, — 1. 
For each e > define 

m ■■= -P { - 



and finally put B^{s) P^{Qd(. The e regularizing problem then becomes 



Minimize 



|£;^(/):= J{{A{X,Df),Df)+TB,{f)}dX\feW''P{0),f\g^ = ^'^, (4.1) 



The existence of minimizers of l|4.1|l is standard. One then proves Lipschitz regularity and 
nondegeneracy for v^, uniform in e. By letting e \ 0, up to a subsequence, v'^ will converge to a 
locally Lipschitz function v that is a minimizer of E^- . We omit the details of the proof of Theorem 



Proposition 4.2. There exists a positive constant 7 > 0, depending only on dimension, A, dD, T and ip 
such that any (possible) minimizer u* of problem i2.13[ satisfies 

D-y := {X e D'^' I dist(X, dD) < 7} C {u* > 0}. 

Proof Let P edDhe fixed and B = Br{Y) c D satisfy 

BndD = {P}. 

By a compactness argument on dD, we can select an r < 5So, where Sq is the universal number 
from l|2.6b , such that the above holds for a.e. P E dD. In view of Theorem I4.2[ there exists a 
minimization, v = v{t) to 

{ Er{f):= I {{A{X,Df),Df)+TX{f>G}]dX\f (,W^'^{bB\B),\ 
Min I \. (4.2) 

[ /Ia5s = 0' AaB = f^^ J 

Here r > is a constant to be chosen later. For future reference, let us label the following the sets 

Q:={X (^D'^ r\hB\v{X)>u{X)} and O -.^ {X D^ r\hB \v{X) > 0}. 

It is important to keep in mind that, from the properties of v, we can ensure that there exist 
constants 0, 5 > depending only on A, r, dD and inf Lp such that 

\OC^D^\> e{T), and dist (P, (90 n D*^)) > <5(r). (4.3) 

We now define the function m : 5B\B ^ IR+ as 

r viX) in iD\B)n5B 

> \ mm{u{X),v{X)} in D^n5B, 

Since m competes with v in the minimization problem l|4.2) l, we have Et{v) < Er{m). Hence, the 
following inequality holds 

/ {A{X, Du),Du)dX - [ {A{X, Dv), Dv)dX > r {iL" (O) - ({u > 0} n O)} (4.4) 
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Our strategy now is to obtain a competing inequality to ||4.4| |. To this end, let us consider the 
function OT: M+ defined as 

m{X) max{?;(X),M(X)}, 

and compare it with u in terms of the minimization problem l|2.13|l . Using the minimality of u, 
we obtain 

Qx(ii''{0) + U'{{u>Q})^L^{{u>Q]r\0)) - eA(i:"({">0})) 

> / T{X,dAu)~T{X,dA^)d7r-^[X). 

JdD 

(4.5) 

From properties 1 and 2 of T, and the Lipschitz continuity of the penalty term g\, we conclude 
from | |4.5| | that there exists a small constant ao = ao{dD, T) such that 



— (£"(0) -£"({?/ >0}nC')) > / {dAU~dATl}dH''-^{X). (4.6) 

"0 JdD 

Applying the Divergence Theorem (see the representation in (13. lb ) and taking into account that 
v{X)hv{X) = a.e., we obtain 

/ {dAU - a^ani dn''~\X) > ^— / (ADu, Du) - {ADv, Dv)dX. (4.7) 

JdD sup ifi Jq 

dD 

As a combination of l|4.4|l , | |4.6| | and l|4.7b we deduce that 
A r 



sup Lp ■ 

dD ao 



£" (O) - £" {{u > 0} n 0) > T (O) - ({w > 0} n 0) 



Thus, if r is chosen big enough, depending only upon dimension. A, dD and (p, there must be 
the case that 

Oc{u* > 0}. 

This together with l|4.3l l ultimately finishes the proof of the Proposition. □ 

In order to advance in our analysis, we need another related free boundary problem: an 
^-obstacle type problem, which again, with the free boimdary technology available, is easy ac- 
complished and therefore we omit the details. 

Theorem 4.3. Let M be a measurable set in . There exists a unique function b, solution to the 
following obstacle-type problem: 

Min |y {A{X,Df ),Df)dX \ f e W^'P{D^) f = tp on dD and f < in . 

Furthermore, supip > b > 0, Lb = in {b > 0} and J bLbdX = 0. 

We now can state and proof our main theorem concerning the existence of an optimal config- 
uration to weak formulations of problem l|2.4| l, namely problems and (*Pa)- 

Theorem 4.4. There exists an optimal configuration Vl\ to problem ( |2.30| > (the penalized problem (*Pa))- 
Furthermore, for a universal modidus of continuity a, the A-potential associated to u\, is a-continuous 
in D^ and WuXWc < K{X, D, ip, F, A). 
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Proof. Before starting the proof, let us explain its strategy. We will initially establish the existence 
of a minimizer to a very weak formulation to problem (|2.101 . Afterwards we "regularize" the 
minimizer we have found via a stabilization phenomenon. Here are the details: Lemma 13.41 
assures, for each < (5 « 1, the existence of a function G satisfying 

3 a (wi) = min^A- 

V{S} 

Furthermore, by noticing that the same computation employed in the proof of Theorem 14.21 ap- 
plies if we restrict ourselves to configurations mV{S), we know that 

D., c {ui > 0}, ys > 0. 

Let B = Br(Xo) be a fixed ball in D'~^ and b be the solution provided by Theorem l4.3l to 

Min I J {A{X, Df), Df)dX \ f e W^'P{D'^) / = <^ on dD and / < in {ui = 0} \ b| . (4.8) 

We also consider [) to be the ^-harmonic function in B that agrees with on B'^ . It is standard 
to verify that 

< It < b < () < sup cp. (4.9) 

As before, (more precisely, as in the proof of Proposition |42| taking into account that J b£,bdX = 
0, we find 

/ r{X,dAu)-T{X,dAb)>ci( [ {A{X,Du{),Du{)dX - ( {A{X,Db),Dh)dX] , {A.IO) 
JdD \Jdc Jd^ J 

for a universal positive constant Ci > 0. However, () competes with b in the obstacle problem 
l lOl l, thus, I I4.10I I becomes 



/ r(X, Dau) - r(X, dAb) >cJ I {A{X, Dui),Dui)dX - [ {A{X, Dl)),D\))dx) . (4.11) 
JdD \Jd^ Jdc J 

For the moment, let us assume p > 2. If we take into accoimt Lemma [3.2i we can enhance the 
estimate by below in l|4.11b as 

/ T{X,dAu)^riX,dAb)>C2( [ \v{ui-i))ixfdx) , (4.12) 

JdD \J_DC- / 

for an appropriate positive but small constant C2. Our next step is to compare and b in terms 
of the fimctional 3a- By doing so, in view of l|4.12|l , we obtain 

A£" {{X e BriXo) I uiiX) = 0}) > C3 (^J^^ |V {ui - f)) (X)]" dX^ , (4.13) 

for another constant C3 > 0, depending on dimension. A, sup tp, and F. If 1 < p < 2, we obtain 

[XU' {{X e BriXo) \ uiiX) ^ 0})Y^\ J^jVui\PdX ' >C3(^I^JV{ui~t)){xfdXy 

(4.14) 

In any case, our conclusion is that if Br{Xo) C D^, then \{X e Br{Xo) \ u{(X) = 0}| =0 and 
consequently, from either | |4.13| | or | |4.14|I , is ^-harmonic there. Of course ) < Ziu^^ ) , pro- 
vided Si < 82. However, from the fact that Lu^ = in we have a much stronger conclusion: 

Z{u{-)=Z{u'^). 
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whenever 5i, 62 < 7- We have proven the existence of a minimizer to C'^^^^^), that is, problem 
ilzTa . 

Our next step is now to prove that O* ;= {u^ > 0} is a minimizer to problem | |2.10| |. For that, 
we have to show 

LuX = OinO*. 

Well, but again it is a standard argument to show from either | |4.13l l or | |4.14|I that u*^ belongs to 
an appropriate De Giorgi's class (recall f) is Holder continuous by elliptic estimates). Therefore, 
there indeed exists a modulus of continuity a {(T{t) = |<|", for some a > 0), such that 

\u(X)-u{Y)\ < CXa{\X -Y\). 

In order to prove that £,u = in {it > 0}, we argue as follows: let Xq G {w > 0} be a generic 
point. By the continuity of u, there exists an tq > such that BrgiXo) C {u > 0}. Therefore, in 
view of | |4.13|I or l|4.14|l , we conclude, as before that 

u = [] in Br„{Xo), 

and the Theorem is finally proven. □ 

5 Existence of an optimal shape to problem (12.411 in low dimen- 
sions 

In this section, upon a technical restriction on the dimension, we will show that the original 
volume constrained problem (|2.4|| admits an optimal configuration. The theory that addresses 
the existence of an optimal design for problem ll2.4b in all dimensions will be developed in section 

m 

Our strategy is based on a limiting analysis on the penalized problem (|2.10l l. For that, we 
initially need a simple lemma. 

Lemma 5.1. There exists a constant C > 0, depending on A, T, D and ip, but independent ofX, such that 
ifu\ is the A-potential associated to an optimal shape ^\for problem (I2.20P , then 



I \Vu\{X)\^'dX <C. 

J DC 



Proof. Let O be your favorite smooth configuration surrounding D that satisfies 

{0\D)= L, 

and let u; be its ^-potential, i.e., the ^-harmonic function in O \ I? taking (p and as boundary 
data on dD and dO respectively. By the minimality property of fi^, we know 

/ T{x,dAul)dn"-\x) < ^xini) 

JdD 

^ ^'f^ (5.1) 
- / T(X,dAU^)dW'-\X) 
J_dD 

= Co, 

where Cq is imiversal, as it depends only on you choice for O. On the other hand, using the 
results and notations of Lemma [3. II we have 



{AiX,Dul),Dul)dX = / uliX)dfi^.{X) 

< sup^.Ai„*(K"\i^) (5.2) 



dD 

sup^- / dAul{S)dn"-\S). 

dD JdD 
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From ellipticity and I I5.2I I, we conclude 

ci J^^ \yuiixwdx < ^„_l^gj^^ J^^ dAui{S)dn--\s\ (5.3) 

where Qi is a positive number that depends on A, (fi and D. Now, for each Y G dD fixed, we 
obtain from II5.3II 



(5) 



(5.4) 



In the last inequality we have used Jensen's Theorem. If we integrate (|5.4| | with respect to Y over 
dD, we reach the following conclusion 

/ ^{y.c^l \Vul(X)\PdX) dH'^-^Y) <C2 f r(X,dAul)dn''-\X), (5.5) 

JdD \ JdC } JdD 

where C2 depends only on dD and the non-linearity T. Finally, if we combine ll5.lt , I l5.5b and 
Il2.5b , we deduce that there must exist a constant C > depending only on A, F, D and ^p, such 
that 

/ \Vu\{X)\PdX < C, (5.6) 

which is precisely the thesis of the Lemma. □ 

Theorem 5.2. Assume the dimension n is less than p. Then there exists an optimal configuration f2* to 
problem ||2.4I >. 



Proof. Because of Lemma ( 15. It , up to a subsequence, we can assume ux converges, as A ^ 00, 
weakly in W^'P{D'-^) to a function u* . Furthermore, since we have assumed n < p, it follows 
by the classical Sobolev Imbedding (see, for instance, |Adams75ll ), that passing to another subse- 
quence if necessary, we can further assume that u\ converges locally imiformly to u* in M" \ D 
and thus, u* is continuous in D'^ . We claim that 

Lm* = in SI* := {X D^ \ w*(X) > 0}. 

Indeed, let Xq G Vl* be an arbitrary point in the set of positivity of u* , say u*{Xq) = > 0. By 
continuity, there exists an ?'o > such that 

u*{X) > ^inB,„(Xo). 
Since u\ converges uniformly to in Br^^ (^o)/ there exists a Aq large enough, such that 

> ^inS,„(Xo), VA>Ao. 

However, we have proven that Lu^ = Oin{u^ > 0}. Therefore, for A large enough, each u\ is 
^-harmonic in Br^ {Xq ) . Thus, as argued in the proof of Lemma 13.41 we in fact conclude u* is 
^-harmonic in its set of positivity and the first claim is proven. 
Notice furthermore that, in view of Proposition l4.2[ 

disi{dD,dQ.*) > 7, 

for some 7 > 0. From inequality l|5.1t , we have, in particular, that 

x{z"{ni\D)-i)+ <Co, 
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for a universal constant Co- Thus, using Fatou's Lemma we see that 

{L" (Q* \ D) - l)+ < liminf(i:"(r2*\i:>)-t) + 
= 0. 

That is, our candidate to an optimal design for problem jlAi , fl*, does satisfy 

in*\D) < L, 

so it competes in problem II2.41 . Our final step is to show that in fact fl* is an optimal configuration 
for problem ||2.4| |. For that, let £ be any competing configuration for problem i2A\ , i.e., \ 
D) < L, and v its ^-potential, that is, v satisfies 

Lu = in C \ D, v = ip on dD, tj = on 9£. 

In particular € competes with u\ in (*Pa)/ problem (|2.10l l: therefore, 

m) ■■= I T{x,dAv{x))drr-\x) 

JdD 

> f T{x,dAu\[x))dn'^-\x) 

JdD 

> / T{x,dAU*{x))dn"-\x) + o{i), 

JdD 

because of the weak lower semicontinuity feature of 5 proven in Lemma 13.41 Finally if we let 
A ^ oo in the above chain of inequalities, the Theorem is proven. □ 



It is worth to point out that Theorem 15.21 gives the existence of an optimal configuration to 
problem l|2.4^ with no regularity whatsoever on the medium. That is, up to this point of the 
project, the operator A has been a general bounded measurable degenerated elliptic map. How- 
ever, it turns out that in order to advance on the study of existence of optimal shapes for problem 
Il2.4b , with no restriction on the dimension, some extra information is needed to perform appro- 
priate perturbations on the optimal designs n\. This will be the contents of the next two sections. 



6 Continuous medium and fine weak geometric properties of 
the free boundary 

In this section we will prove that the free boundary, 50^ enjoys the appropriate weak geometry. 
This feature will allow us to produce geometric-measures perturbations that will ultimately lead 
us to conclude that, if the penalty term A is too large, but still finite, then in fact, obey £" \ 
D) < L. The latter will be carried out in section[7l 

As highlighted in the last paragraph of the previous section, in order to accomplish a deeper 
understanding on the free boundary dft\, we will need to enforce a mild continuity assumption 
on the medium. Thus, hereafter, unless otherwise stated, we shall assume that for some e > 0, 
the map 

eC'(K"\L>), V<eeM". (6.1) 
Mathematically, condition (|6.1l l enables C^ " elliptic estimates for solutions to 

LV- = 

and, at least equally important, it imlocks the Hopf 's maximum principle for ^-harmonic func- 
tions. 
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As for our first theorem in this section, we will obtain optimal regularity for yl-potentials u\ 
associated to optimal configurations of Problem {'^\), that is. Problem l|2.10t . Notice that 
inside the function u*^ satisfies tu^; therefore, it is locally C^ " smooth. However, from the 
Hopf 's maximum principle, reaches the free boundary with a positive slope, thus Vit^ jumps 
from a positive value to zero through the free boundary, dil\. The conclusion is that the optimal 
regularity we can hope for is Lipschitz continuity. This is the contents of the next Theorem. 

Theorem 6.1. Let fi^ be an optimal configuration to Problem (I2.20P and u\ its A potential. Then, 
for a constant C that depends only on A, T, ip and D. 

Proof. We will provide two proofs of this important theorem. The first one follows the glamorous 
approach suggested in [lACSl] . Unfortunately, for non-local problems like ours, the efficiency of 
that method is restricted to the case p >2 and a new and more modern argument is required to 
establish Lipschitz continuity for ^-potential associated to an optimal design with when 1 < p < 
2. The second proof we will present works for all p > 1. 

ist pYQQj: jyig Q^gg p > 2. We shall initially obtain an competing estimate for Inequality | |4.13l l, 
with replaced by u^. Enhancing the notation in the proof of Theorem 14. 4i B ~ Bd{Xo) will 
be a ball centered at a point in fl*^, dist(Xo, 913) ^ c? > dist{Xo,dfl*^) and f) the ^-harmonic 
function in B that agrees with u*^ on dB. For any direction ly G we define 

I i < r < 1 and u^(Xo + dru) = 



if such a set is nonempty; otherwise, we put r,y = 1. For almost every direction i/ the map 
r u\{Xo + drv) is in P 
<1, we can compute. 



r u\{Xo + drv) is in 1]- Thus, taking into account that u\{Xo + dryv) — whenever 



d 

^{Xo + dr^v) = I —{u\ - t)){Xo + drv)dr 



< d-{l~r^y/P'x 



-. i/p (6.2) 

V(()-u^)(Xo+H|fdr , 



where, as usual, p' denotes the conjugate of p, i.e., ^ + y = 1- Now, by the Harnack Inequality, 
we know 

inff)>cif)(Xo), (6.3) 

By 

8 

for a constant ci > that depends only on dimension and A. Here Bv stands for B7^{Xo). Let 
us consider the universal barrier, *B, given by 

diy{AiXo + dX,D^{X)) = in Bi(0) \ St (0) 

OS = 0on9Bi(0) ' (6.4) 
OS = ci in 57(0), 

where ci is the universal constant in ( 16.31 1. By the Hopf's maximum principle, there exists a 
imiversal constant C2 > 0, such that 

nx)>c2(l-\X\). (6.5) 
By the maximum principle and (|6.5l l we can write 

t)iXo + dX) > [)(Xo) • «B(X) > C2()(Xo) • (1 - \X\). (6.6) 
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Combining (|6.2b and l|6.6l l we end up with 



\V{tj-u\)iXo + r,.)\Pdr 



>C3f]^'(Xo)-(l-r,). (6.7) 



Integrating I l6.7b with respect to ly over S" ^, taking into account the definition of r^, we find 



p 



I X{ui=Q}dX <cj |V(f) - u\) {X)\PdX. (6.8) 

J Bd{X„)\Bd,i{Xo) JBa(Xa) 



If we replace, in all of our arguments so far, i?d/4(-'^o) by Ba/A{X), for any X G dBd/2{Xf)), we 
obtain 

fi){Xo)_Y f _X{ui=o}dX <Ci ( \V{i)-u\){X)\PdX, VX e 95^/2 (Xq). 



\ d. J JBi{Xo)\Ba/i(X) ' '■ ' JBi(Xo) 

(6.9) 

Integrating | |6.9l l with respect to X, we prove the following important estimate: 

i^-^X MX eBd{X^)\u\{X)=Q)\<cJ \V{^^ul){X)\PdX. (6.10) 

V " / JBa(Xa) 

Now we argue as follows: let p :~ dist(Xo, 951) and for each < (5 « 1, denote \]s the A- 
harmonic function in Bp+s{Xq) that agrees with u\ on dBpj^siXo). Combining l|4.13|l and (|6.101 
together with standard elliptic estimate, we deduce 

< CX^/P{p + S) + Oil). ^ ' 

Letting 5 \ in 1 I6.IIII we finally conclude 

<(Xo) <cdist {Xo,dn\), 

which clearly implies that is Lipschitz continuous up to the free boundary dil,\ and |1 Vu^ || 00 

2nd p^qqJ The general case. Let us assume, for purpose of contradiction, that there exists a sequence 
of points Xf; e fi^, with 

X.^dni, and ^liM^Z+cx). 
For convenience, we will call := u{Xk) and dk := dist(Xfc, fi^), thus our assumption is that 

^ = 0(1). (6.12) 

For each k, let be a point on dil\ that satisfies 

\Yk-Xk\^dk. 

By replacing Xk by another point Xk, if necessary, because of the weak maximum principle we 
can assume that 

Nk = sup u. (6.13) 

On the other hand, by the Harnack inequality, there exists a universal constant k > 0, for which, 
inf u > KNk- Thus 

sup u > KNk. (6.14) 
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Now for each i < 7 < 1, let f)^ be the L-harmonic function in B^dkO^k) taking boundary data 
equals u*^. By comparing, in terms of the optimal design problem l|2.10t , u*^ and the solution to 
the Obstacle problem in Theorem 14 . 3 1 with M ~ {u\ = 0} \ B^d^ (^fe)/ we deduce, as in the proof 
of Theorem 01 that 

(/ {A(X,Dul),Du\)^^A{X,D\)^),D^^)dx\<\{^dkT. (6.15) 

For each fc > 1, we consider the functions U^, IK^ : B\ (0, 1) given by 

Ml{Z):^^u\[Yu^^duZ) and Jf^(Z) := -^f)(rfc + (6.16) 
From l|614t , we know that 

sup Ml > K. (6.17) 

Bi 

We also know that Oil is the unique minimizer of 

©(«):= / {AiYk+^dkX,Dv),Dv}dX, (6.18) 

JBi 

among functions v e Wq'^{Bi) + 'K], and it satisfies 



^iM{A{Yk+idk,D'Kl)) = in Si 



(6.19) 



A direct computation reveals that 



VUliZ) = ^Vu^(yfe + -idkZ) and similarly ^:KI{Z) = ^t)(yfc + 74-^). (6.20) 



(^(X, ^\^\ • (74)" / (^(i"fe + ldk,DUl), DUl)dX, (6.21) 



Combining ( |6.20t and the Change of Variables Theorem we obtain that 

and the same holds when we replace u\ by and IX^ by "K^. In particular, taking into account 
l(6l2l l, |[6l5l l, we find that 

^ (^(Ffe + 74, i^U^), - (^(n. + 74, i^Oi^), = 0(1) (6.22) 

as fc 00. Furthermore, reasoning as in the proof of Theorem 14.41 we show that the family of 
functions {U^j^^j^ is imiformly continuous in Bi. Now we argue as follows, fix a 7* < 1. From 
the uniform continuity, up to a subsequence, 

Ul IX and J{f ^ M, (6.23) 

imiformly in Bi. Also, Yk Yq. From l l6.19|l and l|6.18t , we obtain that 

dW{A{Yo,D7{)) = inBi 
J{ = U ondBi 

and that 3-C is the unique minimizer of 

■Do(w) := / {A{Yo,Dv),Dv)dX, (6.25) 
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among functions v e VFq^'^(Bi) + However, from (|6.22b , we obtain that 

Do(lX) = Do(^K). (6.26) 
Therefore, U = IK. In particular, U solves the elliptic PDE 

div {A{Yo,DU)) = OinBi. 

However, since 11(0) = 0, by the strong maximum principle, IX = 0, which ultimately contradicts 
1 16.1 7|l and the Theorem is finally proven. □ 

Our next step is to prove that u*^ growths linearly away from dfl*^. Notice that this is the most 
admissible growth rate allowed by the Lipschitz regularity previously proven in Theorem 16.11 
Here is the precise statement: 

Theorem 6.2. There exists a constant c > 0, depending on dimension A, D, T and ip, such that 

A-i/f c ■ dist {Xo, dni) < u\{Xq), 

for any G O^. 

Proof. Let us fix Xq g ^\ near the free boundary and label d dist(Xo, Sfi^). From Theorem 
14.31 there exists a unique solution, 0, to the following obstacle problem 

Min I j {A{X, Df),Df)dX \ f e W^^p{D'^) f = ^ondD and / < in {u\ = 0} U (Xo)| . 

(6.27) 

Recall, in Theorem l4.4[ we proved that u\ is too a minimizer for problem (*P^'''*'*), that is problem 
| |2.13|I and clearly (j) competes with u\ in such a problem; therefore 

/ {v{x,dA<jy)~T{x,dAu\))drr-Hx)>x-^cnd'^, (6.28) 

JdD 

for a dimensional constant c„. Since both u\ and are ^-harmonic in D^, where 7 is the number 
in Proposition 1121 for a constant Ci = Ci(r), we can estimate 

(r(x,9^0)-r(x,a^ul))rfH"-i < I {Oa^ - dAu\) dH'^-^ 

dD Jdn 

< ^ {{A{X, D<P),Dcl,) - {A{X, Du\),Dul)) dX. 

(6.29) 

Here we have used the measure representation provided by Lemma |3A] Now let h satisfy 



L/t, = OinB2^(Xo) \Bd(Xo), = in (Xq), and h = 1 on dB2^i{Xo). 

By the Harnack inequality, there exists a constant C2 > 0, such that 

ul{X) > C2ul{Xo)h{X) in B2^{Xo). (6.30) 

Consider the auxiliary function 

( unn{u\{X),C2u\{X^)h{X)] in ByiX^) 

nl{X) in D^\B.,{X,). 

Notice that g competes with in the obstacle problem, thus, combining l |6.28t , l |6.29t and replac- 
ing by Q, we obtain 

A-^ca < ^ ^ {{^X, Dq),Dq) - {A{X, Du\), Du\)) dX. (6.31) 
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Where set set of integration in the above estimate can be taken to be 

n := e B.,iXo) \ B.,{Xo) I C2uUXo)h{X) < ul{X)} . 

However, in this set, we can estimate 

{AiX,D8),DQ{X)) < A\D5\P 

< \DhiX)\P ■ [c,u\iXo)r 
\c2ul{XoVP 



(6.32) 



In the last inequality we have used the C^ " estimate for h. Finally, a combination of (16.311 1 and 
l|6.32|l leads us to 

X-'^Pcd < uliXo), 

for a constant c = c{n, A, D, T, (p), and the Theorem is proven. □ 

Sometimes it is convenient to express nondegeneracy in any ball centered at a point Xq e O^. 
This is the contents of the next Theorem. 

Theorem 6.3. Let K be a compact set and Xq <E^\r\ K. Then, 

sup u\ > or, 

BAX) 

for some constant c > depending on dimension, K, A, D, T, (p and A. 

Proof. The proof is basically the same of as the proof of Theorem [621 The only difference is that 
u\ is no longer ^-harmonic near a free boundary point Xq, thus we replace the employment of 
Harnack inequality in ( 16.301 1 by: 

v{X) := supit^ • h{X) > Ux on dBr, 

Br 

where h is the ^-harmonic function in \ Br/2 taking boundary data 1 on dB,- and in B,./2- 
We then define the auxiliary function g{X) := mm{u\{X),v{X)}. The proof now follows the 
same path as in the proof of Theorem l6.2l □ 

As usual, optimal regularity, Theorem l6.1l and nondegeneracy, Theorem l6.2l or Theorem l6.3[ as 
you like, allow a deeper understanding on the geometric-measure properties of the free bound- 
ary. In the next Theorem we will show that the free boundary fl*^ has the appropriate weak 
geometry. 

Theorem 6.4. There exists a constant < < 1, depending on dimension, A, D, T, (p, and X^^p, such 
that, 

<;w„r" < {Br{Z) n ni) < (1 - 'i)uj„r^, (6.33) 

for any hall Br{Z) centered at a free boundary point Z G dfly Furthermore, the optimal configuration 
ri^ is a set of locally finite perimeter and for positive constants c, C, depending on A, D, T, ip, and 
there holds 

for any ball Br{Z) centered at a free boundary point. In particular, 7i"^^ (d^x \ 9icd^x) ~ 0. 

Proof. The estimate by below in | |6.33l l, that is, <;a;„r" < iL" {Br{Z) n il*^), is an immediate conse- 
quence of Lipschitz regularity and strong nondegeneracy. 
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Let us focus our effort to prove the uniform density of the zero phase, R" \ rt\. Let us assume, 
for purpose of contradiction, the existence of a sequence of positive real numbers rj with rj \ 

as j ^ oo and 

^ \^ ^ ^=0(1). (6.35) 

We consider then the blow-up sequence : Bi ^R, defined as 

q,{Y):^-ul{Z + r,Y). (6.36) 

Let [}j be the solution to 

diy{AiZ + rjX,D\)j)) = in Bi 



u «R (6-37) 

= on oBi . 

A renormalization of | |4.13|I , when p > 2 or | |4.14|| when 1 < p < 2, imder the assimiption l|6.35b , 
reveals 

/ \V{i),-q,){Y)fdY = 0{l). (6.38) 
By Lipschitz regularity of u\, and C^ " elliptic estimate, up to a subsequence, we may assume 

q/^qo and i), i)o- (6.39) 

uniformly in Bg/n. From I l5.3t [jo satisfies div(^(Z, Dt)o{Y)) = 0, and from Il6.38b so does qo, that 
is, 

div {AiZ, DqoKY)) = in B,^^. (6.40) 

Since q(0) = 0, by the strong maximum principle, we conclude q(0) = in Bi/2- However, this is 
a contraction on the nondegeneracy property guaranteed by Theorem 163] 

We now turn our attention to l|6.34t . The estimate by above, that is id^\ H Br{Z)) < 

(jj.n-1 jg 3 consequence of Lipschitz regularity of u\. In order to prove the estimate by below in 
| |6.34|| , as before, let us assume, for the sake of contradiction, that there exists a sequence rj \ 
such that 

{dn\nBr^{z)) 



7^ . n — 1 



0(1). (6.41) 



With the notation as in l|6.36b , let us define the sequence of nonnegative measures Vj, in -82/3/ as 

j/j div {AiZ + TjX, Dq.j)) dX. (6.42) 

Via a compactness argument, we may assume, modulo passing to a subsequence if necessary, that 
Vj in the sense of measures. However, condition | |6.41| | translates in terms of the measures 
Vj as 

0. (6.43) 

Moreover, by Lipschitz regularity, nondegeneracy and uniform positive density of both phases, 
estimate | |6.33|I , it is not hard to verify that 

Vj ^ z/o := div {A{Z, Dqo)) dX. (6.44) 

Indeed, from I I6.33I I, iL"(9{qo > 0}) = 0, thus in order to justify | |6.44| |, it is enough to attest such 
an identity holds true for balls entirely contained in {qo > 0} and in {qo = 0}. If _B C {qo > 0}, 
then by elliptic estimate, q^ converges to qo in a ^ fashion in B. Thus clearly | |6.44| | is true. 
Now, if B c {qo = 0}, then 

dW{A{Z,D(\n))dx\{B) = 0, 
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so we have to show that i^jiB) ^ as j oo. This is a consequence of nondegeneracy. In fact, 
let B cc B. If there were a subsequence, qj^, for each qj^ ^ in B, then by Theorem 16.31 there 
should exist points Pk^ G B, such that qj^, {Pk^ ) > c > 0. Then, passing to another subsequence, 
Pfcj —> P E B, and since q^^^ converges uniformly to qo, we would reach the conclusion that 
qo{P) > c, which is not possible. In conclusion, if Bk is a nested sequence of balls, with Bk / B, 

then, for some jfc S N, = in S^, for any j > jk- Therefore, i'j{B) ''^^ 0, as desired. 
Having verified | |6.44| |, the observation in ( 16.431 1 tells us that 



div (^(Z,Z?qo)) = 0inB2/ 



3: 



and as argued before, this leads us to a contradiction on the nondegeneracy feature of qo assured 
in Theorem l6.3l □ 

An immediate, yet quite important consequence of Theorem l6.4l is a substantial enhancement 



of Lemma [3.1l for the measure i^u\. 



Theorem 6.5. There exists a Borel function Q\, such that Lu\ = Q\ [dQ\. That is, 



diY iA{X, DuD) (l){X)dX ^ / Qx{S)(PiS)dH"-\S), 

Jdni 

for any £ Cq (R"-\D). Moreover, Q\ bounded away from zero and infinity, that is for a positive constant 
C = C(A, n, A, D, r, </j), there holds 

< < Qa < C < cx). 

As to provide some further insight, allow us to make some loose comments regarding the 
representation Theorem 16.51 The Borel function Qx should be understood as a weak notion for 
the dAu\ along the reduced free boundary ^icd^^A- Indeed, in any peace of 911^, there holds 



h{S) = {A{S,Du\{S)),u{S)), (6.45) 

where v is the unit inward normal vector to dQ.\ at S. However, v{S) = , thus, taking 

into account the scaling feature of A, property (c)(iv), from identity | |6.45l l we reach that 



In a more rigorous way, expression I I6.46I I can be proven to hold in terms of an asymptotic ap- 
proximation, that is, the following is true: 

Theorem 6.6. Let Xo G d^cd^X- Then, for any X <eQ.\ near Xq, we have 

uliX) = BxiXo) {X - Xo, lyiXo))^ + o(|X - Xo|), 

Proof. Indeed, consider a convergent blow-up sequence 

qr{Y) := luliXo + rY) "-^ qo. (6.47) 

Easily, from standard geometric-measures arguments, combined with nondegeneracy and the 
convergence in (|6.47t , we see that 

qo = in {X e M" | {X, iy{Xo)) < O} and {qo > 0} = {X G M" | {X, v{Xo)) > O} . (6.48) 
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Moreover 

div {A (Xo, i?qo) , D(\n) = in {qo > 0}. (6.49) 

Notice that d {qo > 0} is the hyperplane e M" | {X, v{Xq)) = O}: a smooth surface. One veri- 
fies from Theorem 16 ■ 5 1 that 

div {A {Xo, Dqo) , Dqo) = Qa(^o) e R" | {X, iy{Xo)) ^ 0} , (6.50) 

hence, reasoning as before, we reach the following conclusion 

Vqo{Y)-i,{X„) = 0x{Xo), e{{X,i^{Xo)) = 0}. (6.51) 

Recall qo is Lipschitz continuous in the entire M". Let qp be the odd reflation of qo with re- 
spect to the hyperplane {X e K" | {X,v{Xo)) = O}. It is easy to verify that || Vq^ ||ioo(R„) = 
l|Vqo||L-(u,«) < C and that div {A{Xa, Dq*) , Dq*„) = in the whole M". Fr om th e C^-" regu- 
larity of qg, we can employ the beautiful and recent blow-up argument from IKSZI to conclude 
that qQ is an affine function. Thus in view of (|6.51t , we obtain 

qoiX) = 0x{Xo) {X - Xo,iyiXo))+ , 

and the Theorem is proven. □ 

We finish this section by proving the reduced free boundary, dicd^\/ admits a nice "stratifica- 
tion". More precisely, we have 

Theorem 6.7. There exists a collection of hypersurfaces and compact subsets Kj c &j, 

such that 

Furthermore, if X G Kj, the unit outward theoretical normal vector —v{X) to 9rcd^^A zs normal to &j. 

Proof. Let B = Br{XQ) be a generically ball centered at a point of the reduced free boundary. By 
the Lipschitz continuity of u\ and the ellipticity of A, we know there exists a constant L, such 
that 

sup A{X,Dul) < ^. (6.52) 
B 6 

Let T be your favorite nonnegative radially symmetric smooth fimction whose support is Bi . 
Normalize it so that < X < 1; ||X||ii(Bi) = 1. Let be the family of mollification induced 
by X, that is, Ic{X) = e^"2'e(e^^X). Also, select your favorite nonnegative function r] G C^{B), 
satisfying sup?? = L^^ . For sake of notation convenience, let us call V{X) := A{X, Du*^. If v 
denotes the Radon measure Oxn^, we have, for e « 1, 

i^{B) supi / div Vrf^ U' e Co(S;M"), IIV-II < 1 i 



In*. 



> 



I diY{{riV)*I^)dX 

I * div [rjV) dX (6.53) 
div [i^V) dX + 0(1) 

V-V7^dX+ [ QxiS)r]iS)dn''-\S)+Oil). 
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Letting e ^ in l|6.541 and afterwards letting rj L ^ , we conclude there exists a constant 
c(A, A, n, r, ip), such that 

i^iB) > c-W""i (B n dr^d^l) . (6.54) 

hi particular 7Y"~^ [Sred^^A absolutely continuous with respect to Oxn^- Now, arguing as in 
|DeGiorgi55) (see also IIGiusti84l page 54 or IIEG92I page 205) we prove the Theorem. □ 

7 Existence of an optimal configuration for problem (I2.4I > in any 
dimension 

In section |5l upon a restriction on the dimension, we have shown problem l|2.4|l has a minimal 
configuration. The strategy there was to let the penalizing parameter A go to infinity and use 
appropriate estimates that becomes available under the constraint n < p, due to the Sobolev 
Imbedding Theorem. 

The goal of this section is to explore the geometric-measure properties of the free boundary 
dil\, established in the previous section, to settle to existence of an optimal design for problem 
Il2.4b in all dimensions. However, as the readers should expect, the analysis here is rather more 
delicate as we will not be able to pass the limit on the penalty parameter A. Instead, we will show 
that if we adjust the penalty term g\ properly, any optimal configuration, O* = il^, for problem 
II2.IOI1 will obey 

{n*\D) < L. 

Therefore, O* itself will be an optimal design for our primary optimization problem ll2.4b and all 
the regularity features proven to hold for a solution to problem ||2.10| | will automatically extend 
to a solution to problem (12.41 1. 

Before continuing, let us explain our strategy in a bit more technical terms. We will per- 
form a small perturbation on an optimal configuration around a point on the reduced free 
boundary: the portion of 90^ where we can replace classical differential geometry arguments by 
geometric-measures ones. We will not compute the Borel function Q\ of Theorem l6.5l as it is an 
extraordinary hard task: the free boundary condition for problem l|2.10l l is expected to be highly 
nonlocal. Instead, we will show that assuming £"(17^ \D) > l enforces a universal bound to the 
penalty parameter A. 

With the strategy well imderstood, let us establish the first supporting result towards the main 
goal of this section. 

Lemma 7.1. There exists a constant M > 0, depending on dimension, D, ip, V and A, hut independent 
of \, such that 

inf Qx < Af, 
ivhere Qx is the Borel function in Theorem \6^ 

Proof. Indeed, in the lights of Lemma 15.11 there exists a constant C, independent of A, such that 
l^iAllvKi p 1^ C. Thus, from the Trace Theorem for Sobolev functions, we can write 

WulWw^.. ■ {{ul > 0})W > I p>{Z)dH-\Z). 

JdD 

The above estimate combined with the Isoperimetric Inequality assures the existence of a con- 
stant Qi > 0, independent of A, for which the following estimate holds 

n^'-^d^ed^l) > c,. (7.1) 



25 



From l|7.1t and the representation in Theorem l6.5l we have 

(7.2) 



JdD Jdradi^l 



> Ci inf Qa- 



Now, in view of estimate | |5.4|| , for each Y e dD fixed, we establish the following estimate 



dD 



inf Qx 



JdD \ Jdi 



dn''-'{x)< I r(y,4 OauI] dn''~'{x) <C2. 

'dD 



Integrating the above estimate with respect to Y over dD and arguing as before, we conclude the 
proof of the Lemma. □ 

We now pass to describe the mathematical setup for the suitable perturbation technique we 
shall employ on n\ near a point on the reduced free boundary. Initially, select and fix, throughout 
this section, a free boundary point Zq G d^cd^x, such that 

Qx{Zo)<5 inf Qx<Mi, (7.3) 

where Mi depends only on dimension. A, D, T and ip, but it is independent of A. The existence 
of such a point is guaranteed by Lemma [7. II 

Let ifj-.M. M be your favorite nonnegative smooth function whose support equals [0, 1]. 
Normalize it so that 

1p{T)dT ~ I. 

For a fixed positive, but small, real number a, we define the inward perturbation map around Zq 
as 

MX) I ^ - (^T^) ^(^-^^ ^ ^ ^^(^''^ (7.4) 

I X X^BriZo). 

Here, i'{Zo) denotes the theoretical measure outward normal vector at Zq. The idea now is to 
compare n\ with its inward perturbed configuration given by: 

:= {ni) . (7.5) 

For that, let us call Ur the ^-potential associated to fi,-, that is, Ur is the solution to 

hur = in fir \ D 
Ur = (f on dD (7.6) 
Ur = on dflr 

Although it is possible to compare Ur and u directly, it turns out to the more convenient to use 
the auxiliary function, Vr, implicitly by 

VriMX))=uliX). (7.7) 

Notice that {{vr > 0}, w^) is not smtable for our minimization problem | |2.10| |. Also it not efficient 
to compare it with u*^ in terms of the minimization problem | |2.13l l, since = dj(Vr- Our 

strategy is to compare Vr with u\ and with Ur separately and then combine these information 
using Vr as a bridge from Ur and u\. 

The next two Lemmas are from IOT06I , Section 4, though in that paper the computations are 
carried out only for the p-Laplacian operator. Thus we decide to include in this present work 
"economic versions" of their proofs as a courtesy to the readers. 
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Lemma 7.2. With the notation previously set, we have 

{{u > 0}) - {{vr > 0}) = Maar" + o(r"), 
for a universal constant Ah > 0. 

Proof. For sake of notation convenience, we will write u for u\ . For each r > small, we consider 
the r-normalization of u around Zq, Ur'. Bi ^ M, defined as 

Ur{Y) \{Zo + rY). 

Since Zq e 9,.cd^^A' 

Bi n {w, > 0} ^ {y e Si I (y, i.(Zo)) < 0} , (7.8) 

in the sense that the characteristic functions of the above sets in the LHS converge to the charac- 
teristic function of the set in the RHS in the L\^^ (M" ) topology. One easily sees, by the Change of 
Variables Theorem, that 



{Br{Zo) r\ [Vr >())) _ I 



dX 

BAZo)n{v^>a} 
dY 

Bin{v^{Zo+rY}>0} 



I 

J dct{D^r{Zo + rY))dY 



(7.9) 



Bin{ur>o} 



Bin{{Y,iy{Z„))<0} 

It is important to highlight that for any imit vector jy G §"~^, 

J ^'(|y|)^^,;/^rfy^M2, (7.10) 

Bin{(y,i^)<o} 

where A/2 is a constant that depends only on your choice for i/j. Similarly, one finds that 

£"(B,(Zo)n{< >o}) .^0 



dY. (7.11) 

Bin{(y,i/(Zo))<o} 

Combining jV.lOi and (|7.111 , we conclude the Lemma. □ 

Our next Lemma measures the differential on the ^-Dirichlet integral passing from u*^ to v^. 

Lemma 7.3. There exists a constant A/3 > depends on dimension, D, F, ip and ^, but it is independent 
of \ such that 

4- / {{A{X, Dvr), Dvr) - {AiX, Du\),Du\)] dX < aAh + o{a) + 0(1). 



Proof. Again, for sake of notation convenience, we will write u for u*^. Yet for notation conve- 
nience, let us write, for any vector field V , Q{V){X) {A{X, V), V). Applying the Change of 
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Variables Theorem twice and taking into account that Pr maps Br{Xi) diffeomorphically onto 
itself, we can write 

^ ^ e(Dvr)(x)dx = 4r / &{D<!>.r{^;'{x)r^ ■Vu{^;\x)))dx 



BAZo) B,.(Zo 



\ j Q {D^r{Yy^ ■ Vit(r)) X |det {D<^r{Y)) \ dY 

BAZo) 

J e (L'$,.(Zo + rZ)-'^ ■ Vur{Z)) X |det {D'^>r{Zo + rZ)) \ dZ. 



Bin{ur>0} 

(7.12) 

By an explicit computation it is easy to verify that 

D^riZo + rZ)-^ ■ V<(Z) = Vur{Z) + a "^ \^^}\ z, Vur{Z))v{Zo) + o[a). (7.13) 
Furthermore, we can compute explicitly that 

\d.cl{D^r{Zo + rZ))\ = l-a^^^^{Z,v{Z^)). (7.14) 
A straight combination of (|7.12) l, l|7.13t and l|7.14|l , revels that 

^ ^ Q{Dvr){X)-Q{Du){X)dX = -a [ 9 {Dur{Z)) ^ ''jf^^ {Z, v{Zo))dZ + o{a). 



B^(Zo) Sin{ti^>0} 

(7.15) 

It is simple to verify, from the Divergence Theorem, that 



^'{\Z\) 



,Z,u{Za))dZ ^ - J i;{\Z\)dH^-\Z) = I > 0, (7.16) 

Sin{«^>o} Sin{(z,y(Zo))=o} 

with the appropriate integral orientation. Furthermore, by the Lipschitz regularity of u and stan- 
dard geometric-measure arguments we verify that 

{A{Zo +rY,Vur),Vur) Qx{ZQ)i^{Za)xBin{{YMXi))<o}, (7-17) 

in LP{Bi). Thus, letting r ^ in ll7.15l l, and taking into account | |7.16|I and estimate (|7.3t , we 
conclude the proof of the Lemma. □ 

We are ready to prove the existence of an optimal design for problem l|2.4|l in all dimensions. 

Theorem 7.4. There exists a positive number Xq, such that if^\ optimal configuration for problem 
\2.10\ and \D) > l, then necessarily, A < Aq. In particular, there exists an optimal configuration 

for problem (12. 4P and it enjoys all the weak geometric features derived in Section [6] 

Proof. Throughout the proof we fix an optimal configuration Q.\ and assume 

(17* \D)> i. (7.18) 

Initially we recall the variational characterization of the ^-potential u,., namely 



j {A{X, Dur), Dur)dX = min | j {A{X, Dv), Dv)dX \ v^ipondD and u = on arj^j 



(7.19) 
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Now we compare with ilr in terms of the minimization problem (|2.10l l. From the minimality 
feature of the configuration if r is small enough as to iL" {ilr \D) > t, we have 

A{£"(f7* < / r{x,dAUr)-r{x,dAui)dn''-'^ix). (7.20) 

JdD 

As argued before, we have the following estimate 

/ r{X,dAUr) ~T{X,dAul)dn''-\X) < C{dD,r) [ {dAUr - dAul} dW'-'^iX) 

JdD JdD 

< C{dD,r,mfip) I {A{X,Dur),Dur) 
-{AiX, Dul),Du\)dX. 

(7.21) 

Now combining Lemmas[Z2]and[Z]3]with U.19\ , l l7l20|l and \7.11\ , we obtain 

A {A/sar" + o(r")} < C{dD, T, inf ^)r" x [aAh + o(a) + 0(1)] . (7.22) 

If we divide expression Il7.22b by r", let r — > and afterwards divide the result by a and let a \ 0, 
we finally conclude the proof of the Theorem. □ 
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